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(54) Title: SPECTROSCOPIC ANALYSIS OF SAMPLES WITH TURBIDITY AND HIGH ABSORBANCE 



(57) Abstract 

A method of performing analysis of samples having turbidity and/or high absorbance is 
disclosed. The analysis utilizes information from multiple factors, which affect transmittance 
or absorbance.These factors include: absorbance as measured in a non-turbid sample, 
scattering losses due to particles and limitations of the measurement device as measured 
in a non-absorbing sample, additional scattering losses due to sample absorbance, variable 
path length effects due to the sample, variable path length effects due to the measurement 
device, and additional nonlinear effects. By taking into account additional factors that affect 
transmittance or absorbance, a more accurate analysis is achieved. 
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TITLE OF THE INVENTION 
SPECTROSCOPIC ANALYSIS OF SAMPLES WITH TURBIDITY AND HIGH ABSORBANCE 



' 0 CROSS REFERENCE TO RELATED APPLICATIONS 

This application claims priority under 35 U.S.C. § 1 19(e) to provisional patent 
application serial number 60/060,018, filed September 25, 1997; the disclosure of which 
is incorporated herein by reference. 

1 5 STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 

DEVELOPMENT 
Not Applicable 
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BACKGROUND OF THE INVENTION 
The estimation of concentrations of one or more components of a multicomponent 
sample is typically performed by spectroscopic analysis. Generally, the presence of other 
components within a sample makes it difficult to determine the concentration of a single 
component within a sample from a single measurement. Multiple measurements are 
required, and a multivariate technique is used to determine the concentration of a specific 
component from the measurements. Certain assumptions can be made about a sample 
measurement spectrum. Each component of a multicomponent sample has an individual 
absorption spectrum, and the absorbance spectrum of the multicomponent sample can be 
assumed to be the linear combination of the absorbance spectra of each component of the 
multicomponent sample. Further, the absorbance spectrum of each component of the 
multicomponent sample is linearly proportional to the concentration of the component 
within the multicomponent sample. 

The analysis of absorbance in turbid and inhomogeneous samples has proven to be 
difficult, and is further complicated if the sample has high absorbance. A radiation source 
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provides radiation, typically light, through the sample, with a detector detecting the amount 
of radiation passing through the sample. For turbid and/or high absorbance samples, some 
of the provided light is scattered, and may become undetected by the detector. In addition 
detected scattered light may have traversed the sample in a circuitous way, giving rise to 

5 multiple path lengths. All these effects caused by the scattering lead to erroneous readings 
of absorbance factors of the sample being analyzed. Further, some turbid samples may 
consist of relatively large interparticle spaces where no absorbance occurs, thus allowing 
some light to pass through undiminished. This may also result in erroneous readings of 
absorbance factors of the sample. A particular example of this is in the analysis of whole 

1 0 blood for hemoglobin. 

Traditionally, when analyzing whole blood for hemoglobin, the blood must be 
preconditioned to remove turbidity and adjustments made to insure that the measurement 
takes place within the optimal absorbance range of the analysis instrument. The turbidity of 
the blood is reduced or eliminated by disrupting the red blood cells in the sample in order to 

15 release the hemoglobin. The absorbance is adjusted by dilution of the blood or by 
appropriate choice of cell sample cell dimensions. 

Prior attempts to perform analysis of turbid and/or high absorbance samples have 
included utilizing specialized apparatus that collects the scattered light as well as the directly 
transmitted light. Others have utilized an integrating sphere to perform the collection of the 

20 scattered light and the directly transmitted light. A further attempt utilized only pre-selected 
monochromatic wavelengths chosen for a particular application. It would be desirable to 
perform analysis of turbid and/or high absorbance samples without having to precondition 
the sample and without having to utilize specialized light collection devices. 

25 BRIEF SUMMARY OF THE INVENTION 

A method and apparatus for performing an analysis of samples having turbidity 
and/or high absorbance is disclosed. The analysis utilizes information from multiple factors 
that affect transmittance or absorbance of light through the sample. These factors include: 
absorbance as measured in a non-turbid sample, scattering losses due to particles and 
limitations of the measurement device as measured in a non-absorbing sample, additional 
scattering losses due to sample absorbance, variable path length effects due to the sample 
(scatter, bubbles, obstructions and the like), variable path length effects due to the 
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measurement device, and additional nonlinear effects. By taking into account additional 
factors that affect transmittance or absorbance, a more accurate analysis is achieved. 

The present invention, by taking into account non-linear factors as part of the 
detennination of the concentrations of the components of a multicomponent sample is able 
to provide a more accurate estimation of the concentrations of the components of the 
multicomponent sample in the presence of scatter, bubbles, obstructions and the like. 



BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 
The invention will be more fully understood from the following detailed 
description taken in conjunction with the accompanying drawings in which: 

Fig. 1 is a graph of an absorbance spectrum of a normal sample and an obstructed 
sample; 

Fig. 2 is a flow chart of the presently disclosed method for determining analysis of 
15 samples; 

Fig. 3 is a flow chart of the presently disclosed method for determining a residual 
vector for determining analysis of samples; and 

Fig. 4 is a block diagram of a spectrometer of the present invention. 



DETAILED DESCRIPTION OF THE INVENTION 
The analysis of the absorbance spectrum of a sample that is turbid and/or has high 
absorbance can be performed using conventional spectroscopic equipment. The analysis of 
the absorbtion spectrum of a sample is described in terms of absorbance, as opposed to 
transmittance. Absorbance is defined as the ability of a layer of a substance to absorb 
radiation, while transmittance for the purposes of this application is denned as the fraction 
of radiant energy that having entered a layer of absorbing matter reaches its farther 
boundary. The relationship between absorbance and transmittance is defined by the 
equation : 

a = -In t 

Where: a is the absorbance, and 

t is the transmittance 
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The measured spectrum of an ideal sample can be described as: 



a=E 



where: 



a is the column vector denoting the absorbance at particular wavelengths, 
E is a matrix of column vectors, each denoting the spectrum of a 

component or factor at a particular wavelength, and 
c is a column vector describing the particular concentrations of 
components or factors, which may include scatter terms. 

Difficulties occur when attempting to measure samples which are nonlinear, such as 
samples which are inhomogeneous, or when scattering and instrument configuration cause 
the light path through the sample to have multiple path lengths. The averaging of 
transmittance over the multiple path lengths does not result in an averaging of absorbances, 
due to the nonlinear relationship between transmittance and absorbance, as described above. 
Instead the averaging of the transmittance in this case can be characterized by a nonlinear 
function of absorbance, such as cosh(a), which can then be included in the analysis by 
treating the nonlinear term as an additional component. For example the square of the 
absorbance is an excellent approximation for this nonlinear term. For nonlinear terms the 
spectrum becomes : 



where: 
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E„. s [E | /(a) ] 



where y(a) is the nonlinear term(s), and 
where cf {a) is the concentration associated 
with the nonlinear term(s). 



In a specific example, the use of a nonlinear component to detect and correct for 
sample inhomogeneity or instrument effects for lysed blood is shown. In this example, 
sample chamber obstructions produce an inhomogeneous measurement zone that interferes 
with normal optical measurements. The nature of many of these obstructions is to create a 
zone in which blood is excluded. This exclusion results in higher optical transmittance for 
that area. An example of two sequential samples of the same lysed blood sample is shown 
in Fig. 1 . The sample with the higher absorbance gives the normally expected spectrum, as 
indicated by the waveform labeled "Normal". The sample with the obstruction has an 
abnormally low absorbance, as shown in the waveform labeled "Obstructed". The 
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difference between these spectra is shown as the dotted line labeled "Difference". The 
difference is modeled by the square of the normal absorbance multiplied by a negative 
constant labeled "Normal 2 " which almost exactly overlays the "Difference" waveform. 

The fractions and tHb levels predicted by a known method are shown below in 
5 Table 1 . Notice that the tHb level drops and that the fractions change dramatically. 
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TABLE 1 
Previous Method 



tHb Hhb Q2Hb COHh MetHh Scatterl Scatter2 T.v^H 2 - FIT 

14.79 -1.94 102 0.63 -0.4 -0.7 1.63E-04 
12.37 8.17 82.39 3.24 6.2 -1.06 2.23E-03 



1 5 The fractions and tHb levels predicted for the same two samples using the presently 

disclosed method are shown below in Table 2. In this instance the tHb change is less than 
0.3 g/dL and the worst case fraction change is only approximately 2%. It should also be 
noted that the Lysed estimate is a relatively large negative value for the obstructed sample 
as compared with the normal sample. The large negative value is typical for most 

20 obstructions, though some obstructions can cause positive values. 

TABLE 2 
Presently Disclosed Method 



25 




14.49 -0.73 99.91 0.69 0.13 -1.15 0.36 0.00592 7.68E-05 
14.22 -2.22 99.16 2.75 0.31 0.76 -0.58 -0.0428 1.95E-04 

The present method operates by first estimating the fractions (100% 02Hb is 
assumed) then using the square of this spectrum as an additional component in the least 
squares analysis. The square of the new estimated spectrum (based on the resulting 
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fractions obtained in the least squares estimation) is then incorporated and the process is 
then repeated until the fraction values converge. Typically only 2 or 3 iterations are 
required if the first estimate is good. 

The convergence process is illustrated below in Table 3 using the obstructed 
5 sample as an example. 

TABLE 3 

Convergence of the Presently Disclosed Method with an Obstructed Sample 
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14.25 


-1.69 


99.15 


2.25 


0.3 


0.81 


-0.69 


-0.0430 


2.14E-04 


14.22 


-2.2 


99.15 


2.73 
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0.76 


-0.58 


-0.0428 


1.95E-04 


14.22 


-23.2 


99.16 


2.75 


0.31 


0.76 


-0.58 


-0.0428 


1.95E-04 



The definition of the Lysed 2 term is illustrated as follows. When the sample is 
inhomogeneous or when the path of light through a sample is not single valued, such as 
when the light is highly diverging or converging, or when the sample includes bubbles or 
other obstructions, the transmittance is defined as the average of the transmittances of all 

20 light rays. Absorbance, as described earlier, is the negative log of transmittance. For 
simplicity and without loss of generality, consider only two such light rays that traverse the 
sample with different pathlengths. The absorbances of the two rays may be written as a + 8 
and a - 8 respectively, where 28 represents the absorbance difference between the two paths 
taken by the rays. The total absorbance measured is then given by 

25 Abs = -In {e- a+s + e - a -*)l2 [WSGH1] 

This expression can be simplified by removing the common factor e~ a . 



Abs = -In e' a - In (e +s + e-*)/2 



30 



6 



WO 99/15880 PCT/IB98/01492 



Further simplifying the expression and recognizing that the second term contains 
the hyperbolic cosine: 



Abs =a - lnfcosh 8] 



The natural log and hyperbolic cosine can be expanded into a series as follows: 



cosh8=l+8 2 /2!+8 4 /4! + ... 



10 lnl+x = x-x 2 /2-+- x 3 /3-... 



Using these expansions in the expression for absorbance gives: 



15 



Abs = a - 6 2 /2 (1 + 5 2 /12 + 8 4 /360 + ...) + S 4 /8 (1 + 8 2 /12 + 8 4 /360 + ...) 2 + ... . 
For reasonably small values of 8 the above expression can be simplified to : 



Abs = a - 8 2 /2 

20 Using the assumption that the deviation 5 from the average absorbance a is 

typically proportional to a, the above expression becomes: 

Abs = a - k *a 2 
where k is a proportionality constant. 

25 

It should be noted that this expression is in the form of the approximation used to 
model the measured absorbance a. The constant k is determined during the least squares 
analysis. Although the process has a nonlinear term a 2 , the present method uses an iterative 
linear least squares process in which the nonlinear term is treated as an additional linear 
30 term which is fixed for the iterative step. It should also be noted that higher order terms in 
the above expansion may be included in the analysis described herein as needed. 
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The presently disclosed method is shown generally in Fig. 2. The method includes a 
first step 20 in which a matrix of column vectors is established. Each column vector of the 
matrix denotes a spectrum of a component of the sample at particular wavelengths and 
further includes a nonlinear term as part of the matrix of column vectors. The nonlinear 
5 term takes into account the effects of absorbance as measured in a nonturbid sample, 
scattering losses due to particles and limitations of the measurement apparatus, additional 
scattering losses due to absorbance, variable path length effects due to the sample, variable 
pathlength effects due to the measurement apparatus and additional nonlinear effects. 

At the next step 30 a multicomponent column vector denoting absorbance of the 
10 multicomponent sample at particular wavelengths is established. 

At step 40 a concentration column vector denoting estimated concentrations of one 
or more components is determined, such as by dividing the matrix of column vectors by the 
first column vector. The resulting second column vector of estimated concentrations is 
more accurate than prior art methods since the effect of nonlinear factors has been integrated 
15 into the determination of the matrix of column vectors from which the second column 

ft 

vector is determined. 

In certain instances the instrument or sample effects are difficult to determine by a 
priori analysis. For example, variations in the refractive index of the sample, also known as 
anomalous dispersion, cause scatter that varies with the concentration of the sample 
components. This scatter can be characterized by a method that incorporates known 
information about the sample by a hybrid method of analysis. Unknown components of the 
sample or instrument effects can also be considered using this method. 

Hie absorbance spectra of a sample (as measured in a non-turbid sample) can be 
estimated using reference methods to determine the concentrations of the sample 
components and combining this information with known absorbance spectra of the 
components as shown : 

a ref ~ E * C^f 
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where: a ref is the estimated absorbance, 

E is a matrix of column vectors, each denoting the 
spectrum of a component or factor at a particular 
wavelength, and 

c ref is the estimated concentration from reference methods. 
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The estimated absorbance spectrum is combined with other estimated spectra from 
known components with unknown concentrations (e.g. scatter and/or the square of the 
absorbance) and the combined spectra used to approximate the measured absorbance using 
known techniques such as Least Squares Estimation (LSE). 



measured 



where: is the combined spectra matrix, and 

Ctfa, is the vector of estimated concentrations. 

The difference or residual between the estimated spectrum and the measured 
spectrum is the approximation of the unknown component(s). This difference between the 
estimated spectrum and the measured spectrum is used as a vector or vectors in the analysis 
of future samples. Techniques such as Partial Least Squares (PLS) , Principal Components 
Regression (PCR) and Classical Least Squares (CLS) can be used to determine vectors from 
residual spectra. CLS is then used to obtain concentrations of the vectors in a sample. If 
only a single component is known, then the vector for that single component can be 
determined via regression techniques of PLS or Principal Component Analysis (PCA) 
followed by regression testing. 

r exp * * C ui£j - a mtasured 

20 where: is the residual spectrum. 

The analysis using a single residual vector is shown by : 

c new * (E'new * E new )*' * E' new * SL metswti 

25 where: Enew s [E | v res ] v res is (are) the vectors) from the residual 

term(s), and c new is the vector of estimated 
concentrations 



Referring now to Fig. 3 a flow chart of the method 50 for determining a residual 
vector for an unknown sample component or effect is shown. At the first step 60 a 
determination of a first matrix from a first column of vectors of absorbances at particular 
wavelengths is performed. At the next step 70 a second matrix is determined by combining 
the first matrix with an estimated spectra from an unknown concentration or effect. Next, at 
step 80, a third matrix representing a measured absorbance spectrum for the sample is 
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determined. At step 90 a residual vector is determined from the second matrix and the third 
matrix. The residual vector is then utilized as part of the determination of sample 
concentrations in later measured samples. 

A specific example of this is in the measurement of hemoglobin in lysed blood. The 
5 vectors utilized in the analysis are typically the major hemoglobin components of interest 
(e.g. 0 2 Hb, Hhb, COHb, MetHb), hemoglobin components of minor interest (SulfHb and 
CNMetHb), and a term for the scatter or the interference. This analysis by prior methods 
has proven to be inadequate for the accurate measurement of hemoglobin in whole blood. 

A power law of wavelength can approximate losses due to simple particle scatter. 
10 Since the spectral shape of simple scatter as well as its magnitude vary with blood 
components and chemistry, multiple parameters are required for an adequate 
characterization. 

Additional measurement effects must be included for accurate measurements. As 
stated above, hemoglobin is concentrated in red blood cells and the short path lengths 

1 5 required for reasonable optical measurements results in only a small number of cells being 
sampled by the light. This small number leads to results having significant variations in the 
number of cells across the measurements. These variations result in nonlinear absorbance 
variations. As described above, these variations can be approximated by the square of the 
absorbance. The square of the lysed blood absorbance results in better estimates than the 

20 square of the measured absorbance. 

Since hemoglobin is a strong absorber, its absorbance causes variations in the 
refractive index of the red cells. This refractive index variation results in additional scatter 
that is not well characterized by simple particle scatter. Since each hemoglobin component 
has a unique spectrum, the scatter caused by the absorbance depends on the hemoglobin 
25 composition of the whole blood. In the presently disclosed method, this total hemoglobin 
scatter is treated as a linear combination of the hemoglobin scatter from the different 
hemoglobin components. 

In an additional example, an analysis of hemoglobin in whole blood is shown. 
The first step is to determine the vectors that represent hemoglobin scatter. 
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r exp * Etfjj . Ctfjj - a measured 
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Reference values for the samples are then used to determine the individual residual 
vectors by the following: 



Referring now to Fig. 4 a measurement apparatus 100 for performing spectroscopic analysis 
for detemiining estimated concentrations of one or more components of a multicomponent 
sample is shown. The measurement apparatus 1 00 comprises an excitation source 1 1 0 that 
provides an excitation signal to a sample chamber 120. Sample chamber 120 includes a 
multicomponent sample therein. A detector 130 receives the excitation signal that may b 
modified by the sample with in the sample chamber 120. Processing element 140 is 
coupled to, or may be integrated with, detector 130 and determines the concentrations of 
components of the multicomponent sample. This deteimination includes the method 
described above in that the nonlinear effects are taken into account as part of the 
1 5 determination and thus provides amore accurate estimation of the component concentration 
of a multicomponent sample. 

As shown above the presently disclosed method and apparatus provides for more 
accurate analysis of samples having turbidity and/or high absorbance by taking into account 
the effect of nonlinear parameters on the measurement. The nonlinear term is incorporated 
20 into the multivariate analysis resulting in increased accuracy of he analysis. 

Having described preferred embodiments of the invention it will now become 
apparent to those of ordinary skill in the art that various modifications could be made to the 
presently disclosed invention. It is therefore submitted that the invention should not be 
limited to the described embodiments, but only by the spirit and scope of the appended 
25 claims. 
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CLAIMS 

We claim: 



1 1. 

2 



A method for determining estimations of the concentration of one or more 
components of a multicomponent sample using spectroscopic analysis comprising the steps 



3 of: 

4 establishing a matrix of column vectors, each column vector denoting an absorbance 

5 
6 
7 



spectrum of a respective component of said multicomponent sample at particular 
wavelengths and including a nonlinear term as part of each column vector of said matrix; 

establishing a multicomponent column vector denoting an absorbance spectrum of 
8 said multicomponent sample at particular wavelengths; and 

determining a concentration column vector denoting estimated concentrations of one 
or more components of said multicomponent sample based on said matrix and said 
1 1 multicomponent column vector. 

1 2. The method of claim 1 wherein said step of detennining a concentration column 

2 vector comprises dividing said matrix of column vectors by said multicomponent column 

3 vector. 

13. The method of claim 1 wherein said non-linear term is selected from the group 

2 consisting of absorbance as measured in a non-turbid sample, scatter due to sample 

3 absorbance, scatter losses as measured in a non- absorbing sample, variable path length 

4 effects due to said sample, and variable path length effects due to a measurement device. 

1 4. The method of claim 1 wherein said nonlinear term comprises the square of the 

2 absorbance. 



5 . The method of claim 3 wherein said nonlinear term comprises cosh(a), where a 
the absorbance of the respective component of said multicomponent sample. 



is 



1 6. The method of claim 3 wherein said variable path length effects due to said sampl, 

2 comprises variable path length effects due to the presence of a bubble within said sample. 

12 
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1 7. The method of claim 3 wherein said variable path length effects due to said sample 

2 comprises variable path length effects due to the presence of an obstruction within said 

3 sample. 



1 8. 



3 
4 



7 
8 



10 



The method of claim 1 wherein said multicomponent sample comprises lysed blood. 



1 9. The method of claim 8 wherein said nonlinear term comprises variable path length 

2 effects due to said sample including a clot. 

1 10. The method of claim 8 wherein said nonlinear teim comprises the square of the 

2 lysed blood absorbance. 



1 11. 

2 



A method of determining a residual vector for determining estimations of 
concentrations of one or more individual components of a multicomponent sample 
comprising the steps of: 

determining a first matrix from first column vectors of absorbances at particular 

5 wavelengths; 

6 determining a second matrix by combining said first matrix with an estimated 

spectra from an unknown concentration of a component of the multicomponent sample; 

establishing a third matrix representing a measured absorbance spectrum for the 
9 sample; and 



determining a residual vector from said second matrix and said third matrix. 



" 12. The method of claim 11 wherein said step of determining a residual vector 

2 comprises dividing said second matrix by said third matrix. 

113. A measurement device for determining estimations of concentrations of one or 

2 more individual components of a multicomponent sample, the device comprising: 

3 an excitation signal source for providing an excitation signal; 

4 a sample chamber for maintaining a sample therein, said sample chamber in 

5 communication with said excitation signal source; 



13 



WO 99/15880 PCT/IB98/01492 



6 a detector operative to detect said excitation signal after said signal has been passed 

7 through said sample chamber, and 

8 a processing element in communication with said detector, said processing element 

9 operative to provide concentration estimations from data received from said detector, said 

1 0 concentration estimations reflecting effects of at least one nonlinear term. 

1 14. The measurement device of claim 13 wherein said nonlinear term is selected from 

2 the group consisting of absorbance as measured in a nonturbid sample, scatter due to sample 

3 absorbance, scatter losses as measured in a nonabsorbing sample, variable path length 

4 effects due to said sample, and variable path length effects due to a measurement device. 

1 15. The measurement device of claim 13 wherein said nonlinear term comprises the 

2 square of the absorbance. 

1 16. The measurement device of claim 13 wherein said multicomponent sample 

2 comprises lysed blood. 

1 17. The measurement device of claim 16 wherein said nonlinear term comprises the 

2 square of the lysed blood absorbance. 

1 18. The measurement device of claim 1 3 wherein said nonlinear term comprises higher 

2 order terms in the expansion of the absorbance of said multicomponent sample. 
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